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Color quadrupoles have been found to be important in the proper description of observables sen- 
sitive to the small- j: regime in nuclei as well as in the operator definition of the Weizsacker-Williams 
gluon distribution. In this paper, we derive the small-a; evolution equation of the quadrupole and the 
Weizsacker-Williams gluon distribution without taking large Nc limit and study the properties of 
the equation in both dilute and saturation regime. We find that the quadrupole evolution follows the 
BFKL evolution in the dilute regime and then saturates in the dense region due to nonlinear terms. 
■ This leads us to conclude that the Weizsacker-Williams gluon distribution should obey the same 

' geometrical behavior as the dipole gluon distribution as found in the inclusive DIS measurement. 
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^ • I. INTRODUCTION 

(N 

Recent studies of two-particle production processes [J-Q where a dilute probe scatters from a dense target have 
. shown that in order to understand saturation phenomena and the corresponding small-x dynamics it is necessary to 
Q^' go beyond color dipoles and compute the appropriate evolution equations for color quadrupoles as well. In particular, 
the measurement of two-particle correlations in d-f Au collisions [y, 0| has attracted a fair amount of attention due to 
the observed suppression of the far-side peak as the rapidity of the observed particles is increased. So far, the only 
' satisfactory explanation for this phenomenon, showing suppression for forward rapidities and not in the central rapidity 
region, is the calculation presented in which is performed in a Color Glass Condensate (CGC) framework. Even 
though the results of [1] are shown to agree with the current data, it has been argued d, [ll] that the approximation 
^ I used to implement the small-a; evolution of the relevant correlators in terms of only color dipoles is not correct even 
psj ' in the large- A'^c limit. Based on the study of two-particle production processes in DIS and pA collisions [1, |^ it is 
; now believed that for scattering processes involving a dilute probe interacting with a dense target can be described 
C ' in terms of only color dipoles and color quadrupoles in the large- Ac limit. 

' [ , Color quadrupoles have also been shown to be relevant when unintegrated gluon distributions are considered in 
00 the small-a; regime. In small-a; physics, two different unintegrated gluon distributions [8, 9] (also called transverse 
momentum dependent gluon distributions), namely the Weizsacker-Williams gluon distribution xG^^^ [ll| and 
^""^ ■ the dipole gluon distribution xG'^'^\ have been found and widely used in the literature. The Weizsacker-Williams 
" , [ ' gluon distribution, known as the conventional gluon distribution, gives the Fock space number density of gluons 
^ , inside dense hadrons in light-cone gauge. The dipole gluon distribution, defined via the Fourier transform of the color 
' dipole amplitude, has been studied thoroughly since it appears in many physical processes [l^, [13 • This dipole gluon 
, distribution can be probed directly in photon-jet correlations measurement in pA collisions. Recently, studies on 
the Weizsacker-Williams gluon distribution indicate that it can be directly measured in DIS dijet production and its 
operator definition is related to quadrupoles instead of normal color dipoles. Other more complicated dijet processes 
in pA collisions (e.g., qg or gg dijets) involve both of these gluon distributions through convolution in transverse 
momentum space. 

In terms of the operator definitions of the Weizsacker-Williams gluon distribution xG^^-' and the dipole gluon 
distribution a;G'^',the evolution of xG*^^-' appears to be given by the evolution of the dipole scattering amplitude 
which obeys the Balitsky-Kovchegov equation [3 EBl; while xG^^-' is governed by the evolution of both dipole and 
quadrupole scattering amplitudes. For more complicated dijet processes in pA collisions one needs to use both 
the dipole and the quadrupole evolution equations since the cross sections of those processes involve both gluon 
distributions and their combinations. Knowing both dipole and quadrupole forward scattering amplitudes on a target 
allows one to understand their small-x behavior and thus compute physical observables in high energy scatterings in 
large Nc limit. This is necessary to describe dijet processes in both DIS and pA collisions systematically. 

The objective of this paper is to derive the small-x evolution equation of the quadrupole and the Weizsacker- 
Williams gluon distribution without taking the large Nc limit, and study their small-a; behavior. As far as the 
quadrupole amplitude Q is concerned, it seems that qualitatively it is understood quite well. In the leading twist 
limit when the density is low, it obeys a BFKL type equation. When all the sizes of the quadrupole are comparable 
to the scale at which unitarity sets in, namely Axi ~ Qs\y) (QsO^) is the so-called saturation momentum, which 
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depends on the rapidity Y, and whose inverse gives the order of magnitude of the size above which color-neutral 
objects are absorbed.), the quadrupole amplitude evolves towards the stable fixed point Q — and the corresponding 
Weizsacker- Williams gluon distribution xG'^^^ starts to saturate. 

The rest of the paper is organized as follows. In Sec. II, we carry out the derivations for the evolution equation 
of the quadrupole amplitudes and the Weizsacker- Williams gluon distribution. We investigate the weak interaction 
limit of these equations in Sec. III. Sec. IV is devoted to the discussion on the unitarity limit and fixed points of the 
quadrupole evolution equation. The summary and further discussions are given in Sec. V. 



II. THE EVOLUTION EQUATION OF THE QUADRUPOLE 

In general, the JIMWLK evolution [Hi of an arbitrary operator is given by 



with 

^t^\j iu^,v±,zi^) [l + V^K - V^V. - V-^k] , (2) 

where JC {u±, v±,z±_) = (u^-z^yHy^-z^)'^ ^ represents the Wilson line in the adjoint representation. 

If one sets O = (Tr [U {x i)U\y iS)) ^ with f/(xj_) being the Wilson line in the fundamental representation, one 
can easily reproduce the well-known evolution equation for the dipole amplitude, 

— (Tr [U{x)U^y)\}y = -—^ / d z 



dY^ ^ ' ' ^ '^'^ 27r2 ] ^ {x^~ zx_Y{z^-y^f 

X |(Tr {U{x)\j\y)\)^ i- (Tr {U[x)U\z)\ Tr {U{z)U\y)\)^ . (3) 

This is not a closed equation for the dipole amplitude since the last term involves a correlator of four Wilson lines. 
Using a mean field approximation for a large nucleus and in the large- TV^ limit, this correlator can be factorized as 
a product of two dipole amplitudes yielding the well-known Balitsky-Kovchegov equation. The Balitsky-Kovchegov 
equation has a probabilistic interpretation in terms of dipole splittings with the probability of one dipole 

splitting into two new dipoles and {z^,y^) being Pd^dd = {Jilt^)-i'(^l-y^)-i ■ 

Here we would like to derive the evolution equation of quadrupoles by using the JIMWLK Hamiltonian method 
without taking the large Nc limit. By setting O ^ ^ (Tr {U {x\)U'^ {x'-^U {x2)U'^ {x',2S)') y ^^"^ using the above Hamil- 
tonian, we get (see the derivation in Appendix. |X| 

^{Tt{U{x^)U\x'^)U{x2)U\x'^)\)^ 

^ d^z^Kx{x^,x\,X2,x'^-z){TY{U{x{)U\x\)U{x2)U\x'^)\)^ 



i2n) 

+^ J d'z^A{x^,x[,X2,x'^;z)l-{Tr[U\x[)U{x2)]Tr[uHx',)U{x^)])^ 
+ ^ J d'z^B{x,,x[,X2,x',;z)l-{Tr[U{x^)UHx[)]Tv[U{x2)U\x',)])^ 
+^ / d'z±IC2{xv,x[,x',;z)l-{TT[U{x^)UHz)]Tr[U(z)UHx[)U{x2)UHx'^)])^ 
+ ^ J d'zM^'i;^i,^2;z)^^{Tr[U{z)U\x[)]Tr[U(x,)UHz)U{x2)UHx'^)])^ 



asNc 



d^z^IC2{x2;x[,x'^;z)^{TT[U{x2)U^{z)]TT[U{xi)U\x[)U{z)uHx'^)])^ 



(27r)2 J " v-'-i'-.'-.^^ 
^ J d^z^lC2{x'2;Xi,X2;z)^^{Ti[U{z)U\x'2)]Tr[U{x,)U\x[)U{x2)U\z)])^, (4) 
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where 

' ^' ' ^' (x\ — z)"^ [z ~ xyY' (x'2 — z)'^{z — X'l)'^ {x[ ~ z)'^{z ~ X2)'^ {x'2 ~ z)'^{z — XlY'' 

Mr. t' to r'-z^^ « ^ ^2? , {xi-X2f K ^ ^2? (4 " ^1)' 

^ 2'^^ [x'^-zfiz-x'^Y {xi-zY{z-X2Y [x'^-zfiz-x^Y [x'2-zf{z-xiY' ^ ' 

„, , , X K - 4)^ , {xi-X2f (4 - 3:^2)^ K ~ 

D^xi,Xi,a;2,X2,zj [x'^^zYiz^x'^Y {xi-zf{z-X2Y [x'^ - zY {z - X2Y [x'^- zY{z-xiY' ^ ' 

i^2KXi,x^,X2,z, - ^xr~zY{z-x\Y {xi - zY{z - x'^Y {x\- zYiz-x'^Y' ^ ' 

We note that the kernels satisfy the sum rule 
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lCi= A + 13 + Y,^2, (9) 

where the sum in Eq. ([9]) goes over the missing coordinate in K.2 which in Eq. (|8|) is X2- AH the above coordinate 
variables are implicitly assumed to be two-dimensional. Eq. Q El suffers the same problem as Eq. Q in the sense 
that it is not a closed equation because the right hand side includes higher-point correlations. The way to deal with 
this difficulty is the same as for the Balitsky-Kovchegov equation assuming that, for a large nucleus, these correlators 
can be factored as products of correlators involving only one trace at a time when the large- A^c limit is taken. The 
resulting equation is equivalent to the quadrupole evolution equation found in Ref. Q which was derived considering 
only the leading contributions from the beginning. In other words, the full evolution equation has no terms 
which are explicitly suppressed by powers of l/iVj, and taking the large- A'c limit only has the effect of allowing the 
aforementioned factorization of the correlators. This was also the case for the Balitsky-Kovchegov equation [151 as 
can be seen from Eq. 

Let us also comment on the physical interpretation of those seven terms on the right hand side of Eq. ([4]) . The first 
term stands for the virtual correction which subtracts probability from the original quadrupole. The second and third 
term represent the splitting of the quadrupole into two new dipoles while the last four terms indicate the splitting of 
the original quadrupole into a quadrupole and a new dipole. The kernels satisfy the relation Eq. ([S]) as a result of 
conservation of probability. 

Now, let us turn our attention to the unintegrated gluon distributions mentioned in the introduction. First consider 
the dipole gluon distribution whose operator definition is related to the dipole amplitude, namely, the two point 
functions of Wilson lines ^Tr (C/(x^)[/t(y^))^ as follows 

= / |^^"''"'^'""'"^]^(TrC/(x^)[/t(y^))^. (10) 

Its evolution is therefore provided by the Balitsky-Kovchegov equation. 

The operator definition of the Weizsacker- Williams gluon distribution is given by a slightly different operator 
which can be obtained from the quadrupole correlator. According to Ref. [1, [aH^, the Weizsacker- Williams gluon 
distribution can be written as 

xG^^\x,k^)^-- [ ^f^e-^^-<^-^'^i:v[d,U{v)]U\v')[d.U{v')]UHv))y. (11) 

The evolution equatior0 for the correlator (TT[diU{v)]W{v')[diU{v')]U^{v)')y can be obtained from Eq. (jH) by 
differentiating with respect to Xi and X2, and then setting xi = X2 = v and X2 = x[ = v' . Then the resulting 



^ These equations have also been derived using the JIMWLK formalism by D. Triantafyllopoulos and by J. Jalilian-Marian (private 
communications) . 

^ Inspired by Ref. [20{. we find that the small-a; evolution equation of the Weizsacker- Williams linearly polarized gluon distributions is 
given by a slightly different correlator (Tr {diU{v)] {v') [djU{v'y\ {v)')^. We will leave the detailed derivation and phenomenological 
study of the linearly polarized gluon distributions to a future work [2l| . 
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FIG. 1. Illustration of two-step evolution which generates the quadruple amplitude. The dotted lines indicate the moments of 
the interaction with the target nucleus and the dashed line represents the cut. The two dipoles correspond to the two internal 
color lines, and are characterized by the coordinates {zi,v) and (v',22) respectively at the time of the interaction. The single 
external color line interacts as a quadrupole defined by the coordinates {v, z\, Z2,v'). 



evolution equation becomes 



d 



— {Tr[d,U{v)]U\v')[dMv')] U\v))y 



' 27r2 

OsNc 

' 27r2 



{v-v'y 



[v — z)'^{z — w')2 



(Tr [d.,U{v)] uHv') m{v')] U\v))^ 
{v-v')i {v - Z)i' 



(v — — z) 



X {(Tr [U{v)UHv') [d,U{v')] C/t(z)] Tr [U{z)UHv)])^ - (Tr [U{z)U^{v') [d,U{v')] UHv)] Tr [U (v)UHz)]) ^} 



27r2 



cfz. 



{v - v') 



Nc {v- z f{z-v') 



{v' — f )2 (v' — z) 



X {(Tr [[d.,U{v)] U\z)U{v')U\v)\ Tr [U {z)UHv')\)y - (Tr [[d,U{v)] U\v')U(z)U\v)\ Tr [U {v')U\z)\)^] 



a,Nr. 



(fz I 



Nc {v - z)^{z -v')^ 



1 - 



2{{v~z)-{z-v')y 

{v — z)^(z — W')2 



X {(Tr [U{v')U^{z)] Tr [U(z)U^v')])y + (Tr [Uiv)UHz)] Tr [U {z)U^ (v)]) , 
- (Tr [U{v')U\v)] Tr [U{v)U\v')] )^ - N^} . 



(12) 



Among these four terms in Eq. (|12p. the second and third terms are quite troublesome since they introduce new 
correlators involving three coordinates. The first term can be understood as the virtual correction as analogous to 
the first term in the Balitsky-Kovchegov equation shown in Eq. 1^. The last term is in agreement with the results 
obtained from the one-loop calculation in Ref. 22 1. The calculation in Ref. [l^ employs an effective theory which 
allows one to probe the Weizsacker- Williams gluon distribution by using the current j = — 1-^°^^^°^^. One can also 
show that this calculation is equivalent to our above derivation for Eq. (1121) , and one can obtain the first three terms 
as well if all the graphs are included without angular average. 

In addition, we find that one will inevitably run into the evolution of quadrupoles irrespective of the initial conditions 
when the Weizsacker- Williams gluon distribution appears in the process before small-x evolution is included. As far as 
the Weizsacker- Williams gluon distribution is concerned, one finds that after two steps of evolution, the contribution 
of quadrupoles appears as the following S'-matrix amplitude 



^Tr [U{zi)U\z2)U{v')uHv)] ^Tr [U{v)U\z^)] ^-Tr [U{z2)U\v')] 



(13) 



Starting from the correlator (Tr [diU{v)] U^{v') [diU{v')] U''{v))y, the first step of evolution is given by Eq. ([T^ which 
generates terms like (Tr [U{zi)U^ {v') [diU{v')] U^v)] Tr [L''(i')C/^(zi)] ).^. One can further evolve such object and find 
that the second step of evolution yields a combination of a quadrupole plus two dipoles as in Eq. In terms of the 

dipole picture and 't Hooft's double line notation, we illustrate the above terms in Fig. [T] where the double lines at v 
and v' have the highest longitudinal momentum, the double line at zi has the next highest longitudinal momentum and 
the double line at Z2 has the smallest longitudinal momentum. Fig. [T] and its radiative corrections are characterized 
by the feature that the double line in the middle does not directly connect to the virtual photons. By squaring the 
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production amplitude, one gets the basic color and spatial structure of a quadrupole j^Tr [U{zl)U^Z2)U{v')U'' (v)^ 
together with two dipoles j^Tr [?/(«)[/''' (zi)] and j^Tr [U{z2)U^v')] before further evolution. This shows that the 
Weizsacker- Williams distribution does not have a closed evolution equation on its own and, despite its apparently 
simpler structure in terms of only two coordinates, the full quadrupole evolution is needed to include small-x effects. 
In terms of the DIS dijet process considered in 0,0], the correlation limit taken to avoid the quadrupole and be able 
to express the cross section in terms of the Weizsacker- Williams distribution does not help when small-x evolution is 
considered. The same argument and conclusion also holds for the single gluon production in the gedanken DIS process 
using the current j = —\F^uF^y However, for inclusive and semi-inclusive DIS, one can show that the evolution 
always involves only dipole amplitudes and this allows one to write a closed equation, namely the Balitsky-Kovchegov 
equation to describe these processes. 



III. THE WEAK INTERACTION LIMIT 



In order to study the evolution of color quadrupoles in the dilute regime, where the evolution will be described 
entirely by the BFKL evolution equation, let us consider the leading twist approximation to the dipole and quadrupole 
amplitudes and the corresponding evolution equations. In this leading twist approximation, the dipole amplitude 
S{xi^X2)y — ]^(TrC/(a;i)C/^(a;2))i' can be written as S{xi,X2)y = 1 ^ ^^{xi, X2)y with T satisfying the dipole 
form of the BFKL equation 



^r(a;i,a;2)Y = ^^^^ / (fz- -[r{xi, z)y + T{z,X2)y ~T{xi,X2)y] ■ (14) 

aV ZTT^ J (xi — zY[x2 — zY 



For the quadrupole amplitud^ x'^, a;2, a;2)y = ■^(^t:U{xi)U\x'i)U{x'2)W {x2))y , we expand to second order 

in gA with A being the background gauge field. A simple calculation shows that this expansion can be written in 
terms of the linear dipole interaction F (see Eq. (B21) in Ref. ^5,] jf| 



Q{xi,x\,x'2.,X2)y = 1 Y [r(a;i,X2)y -I- T{x[,x'2)y + T(xi,x'i)y -l-r(a;2,X2)i' - F(a;i,a;2)y - T{x\,X2)y] ■ (15) 



For the full quadrupole amplitude, all possible coordinate pairings have to be considered but it is important to remem- 
ber that for the expression of the cross section for the dijet production in DIS, and therefore also for the derivation 
of the Weizsacker- Williams gluon distribution, the quadrupole appears through the combination (5(a;i, x'j^, Xj, 0:2) — 
S{xi^ X2) — S{x'i, X2) + 1. In the linear regime, the contributions from interactions hooking both gluons to the {xi, X2) 
dipole or both gluons to the {x'i,x'2) dipole cancel out. 

Eq. (jl4p above can be obtained directly from the BK equation by expressing the dipole amplitude in terms of F 
and neglecting the non-linear terms in F's. The same approach can be easily followed for the quadrupole evolution 
equation. First, let us rewrite the quadrupole evolution equation in Eq. Q by grouping the different terms with the 



^ Note that we have interchanged X2 and in the definition of the quadrupole in order to match the physical observable calculated in 
the DIS dijet process. 

* The BFKL limit has also been found by D. Triantafyllopoulos (private communication). 
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same dipole kernel 

— Q(xi,Xi,X2,X2) = ^(2:;;^ J d Zj_ 

{Xi - X2f 



(X\ — Z)^{X2 — zY 

(x\ — x'2)'^ 



[Q[xi,x'i,X2tX2) - Q[z,x'i,x'2tX2)S{xi,z) - Q[xi , x'^, x'2, z)S{x2t z) + S{xi , X2)S{x'i, x'2)\ 



{x'l — z)^{x'2 — z)^ 

{xi~x'if 

{xi — zY{x'i — z)2 

{^2 " 4)^ 

{X2 — z)'^{x'2 — zY 
{Xl - X'2f 

{X\ — zY{x'2 — zY 
{x'l - X2f 

{x'l - ZY{X2 - zY 



{Q{xi,x'i,x'2tX2) - Q{xi, Z , x'2, X2)S{x'i, z) - Q{xi, x'l, Z , x'2)S{x'2t z) + S {xi , X2)S{x'i, x'2)\ 

[Q{xi,x'i,x'2,X2) - Q{z,x'i,x'2,X2)S{xi, z) - Q{xi, z, x'2,X2)S{x'i,z) + S{xi,x'i)S{x2,x'2)] 
[Q{xi,x'i,x'2,X2) ~ Q{xi,x'i,z,X2)S{x'2, z) - Q{xi,x'i,x'2,z)S{x2,z) + S{xi,x'i)S{x2,x'2)] 
[S{xi,X2)S{x'i,x'2) + S{xi,x'i)S{x2,x'2) - Q{xi , z , x'2, X2) S {x'l , z) - Q{xi,x'i,x'2,z)S{x2, z)] 
[S{xi,X2)S{x'i,x'2) + S{xi,x'i)S{x2,x'2) - Q{z , x'l, x'2, X2) S {xi , z) - Q{xi,x'i,z,X2)S{x'2,z)] 



(16) 



where we have assumed that, in a large nucleus, one can factorize (0102)^ into (Ci)y (02}y Then we use the 
expressions quoted above for Q and S in terms of F and keep only the linear terms. It is easy to see that each of the 
lines of Eq. becomes a BFKL equation of the form given in Eq. ([Ti| for each pair of coordinates. This is in 

complete agreement with what one would expect from just evolving Eq. (jl5p with BFKL in the linear regime. BFKL 
evolution generates the exponential growth of the quadrupole T-matrix {Tq = 1 — Q) in terms of the rapidity Y. The 
saturation scale Qs{Y) is usually used to characterize the scale at which T approaches 1. If one fixes the relative size 
of all the spatial coordinates of the quadrupole, one expects that quadrupole evolution in the dilute regime should 
obey one single BFKL equation and thus exhibit the same geometrical scaling behavior as for dipoles [23- 25] in terms 
of the variable r^Ql{Y) with r being the typical scale of the quadrupole. 

Similarly, in the dilute regime, the correlator involved in the calculation of the Weizsacker- Williams distribution 
can be written in terms of F's as: 

{Tv[d.,U{v)] U^v') m{v')] U\v))^ = ^0^^0^,r{v,v'). (17) 

Naturally, expanding the correlators in Eq. (jl2l) and keeping only terms linear in F or its derivatives leads to the 
same result as differentiating twice Eq. (|14l) . 

In the dilute regime where the gluon density is low, we know that the Weizsacker- Williams gluon distribution xG*-^' 
and the dipole gluon distribution xG'^'^^ both reduce to the same leading twist result. To see this connection explicitly 
recall that xG^-^"^ is given by the Fourier transform of the correlator in (fT7|) and therefore the derivatives can be 
replaced by fc^ after integrating by parts. 

IV. THE UNITARITY LIMIT AND FIXED POINTS 

Let us begin with the discussion on the dipole amplitude. Schematically, the Balitsky-Kovchegov equation can be 
cast into 

■^S = J Pd^dd{SS-S), (18) 

with S being the dipole amplitude and Pd^dd > being the probability of a dipole splitting into two new dipoles per 
unit rapidity. 

With Eq. (|9]) and the assumption of factorization of operators, the evolution equation of quadrupoles can be 
schematically written as 



■^Q = J Pg^gd {QS -Q) + J Pc^dd {SS - Q) , 



(19) 



where Q represents the generic quadrupole amplitude, Pq^qd is the probability that the original quadrupole splits 
into a quadrupole and a dipole due to gluon emissions and Pq^dd stands for the probability of the original quadrupole 
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splitting into two dipoles. The probabihstic interpretation of the quadrupole evolution equation is of course only 
heuristic since Pq^qd and Pq^dd are not positive-definite. However, we find that Pq^qd + Pq^dd — ^^j^p-^^i is 
always positive (this equality is the identity in Eq. (jO])). It seems to make sense to speak of the probability that the 
original quadrupole decays to qd and dd inclusively, but not that it decays to qd or dd exclusively. 

S = {) and S = \ are two fixed points of Eq. (|18l) . and Q = and Q = 1 are fixed points of Eq. with 
S — Q ~ Q a. stable fixed point and S* = Q = 1 an unstable fixed point. Let us first write S' = -t- (55, Q — Q + 6Q 
where SS and SQ are small perturbations, then Eqs. (|T8)) and (IT9t become -^SS = — J Pd^ddSS and -^5Q — 
— J {Pq^qd + Pq^dd) ■ Siucc Pd^dd > aud Pq^qd + Pq^dd oc /Ci > 0, the cvolutiou cquatlous will drive SS and 
SQ to zero, thus S — Q — is a. stable fixed point. Now let us take S = 1 — SS and Q — I ~ SQ, then Eq. ^TE\\ 
becomes -^SS = J Pd^ddix, y, z) {SS{x, z) + SS{y, z) — SS{x, y)) which is exactly the BFKL equation. Similarly, as 
shown in Sec. HI, SQ obeys the BFKL equation in the weak interaction limit as well. It is well-known that BFKL 
evolution will drive SS and SQ away from zero in terms of exponential growth. Therefore, the fixed point S = Q = 1 
is unstable. We do not know if these fixed points are the only fixed points for the quadrupole evolution, but it is 
certainly true for the Balitsky-Kovchegov equation. As a result, we expect that the evolution will always drive the 
dipole and quadrupole amplitudes from the dilute limit towards the stable fixed point S = Q — Q which leads to 
unitarity and saturation. 

Finally, we note that there is a traveling wave picture [1^, [l^l for the evolution of Q, or of Tq = 1 — Q, exactly 
analogous to that for 5, or T = 1 — 5. The velocities of the traveling waves for Tq and T are identical, since the 
velocity is determined by BFKL evolution. If one scales all coordinates in Tq uniformly then the shapes of the 
traveling wave front of Tq and T are identical except near the top of the fronts, that is where nonlinear terms in the 
evolution become important. 



V. CONCLUSION 

In summary, using the JIMWLK Hamiltonian, we have derived the evolution equation for the quadrupole and the 
Weizsacker- Williams gluon distribution at finite Nc- We find that they follow BFKL evolution in the dilute regime 
and reach the saturation regime as a stable fixed point. Following the discussion in Ref. (27l . [28j . we know that 
BFKL evolution together with a saturation boundary are responsible for the geometrical scaling behavior [23l - [25| of 
the dipole gluon distribution. Since we also observe the same properties for the quadrupole evolution equation, we 
believe that the Weizsacker- Williams gluon distribution should exhibit geometrical scaling behavior as well, although 
its evolution equation is much more complicated. It seems that quadrupoles evolve essentially the way dipoles do, but 
their evolution is more difficult to evaluate due to the complicated structure of the evolution equation. The difference 
of their small-x evolution behavior lies in the transition region between the scaling regime and the saturation regime. 
Likely, numerical studies of quadrupole evolution will be necessary to understand the details of their different behavior. 
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Appendix A: Derivation of the evolution equation for the quadrupole amplitude 

Here we present some essential details of the derivation leading to Eq. Q . We begin with Eq. ^ with the operator 
O set to be ^Tr {U{xi)U'^ {x'i)U{x2)U'^ {x'2)))y- Since there are four different terms inside the square bracket in 
the definition of rj^^^ as shown in Eq. ([2]), we compute the contributions of these four terms separately and provide 
the results as follows. 
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The first term which involves the identity matrix in color space S yields the contribution Ii 



h = 



' 27r2 

asNg 
'(27r)2 



(4 - 



{Tt[U{x,)U\x',)U{x2)U\x'^)])^ 



{x'l — z)'^{z — xiY {x'2 — zY(,z — X2Y 
d'z^A{xi,x\,X2,x'^-z)\^{Tv[u\x'^)U{x2)\TT[u\x'^)U{xi)\) 



^{Tr[U{x,)U\x\)U{x2)U\x'^)])^ 



(Al) 



With the help of the identity T^V^°- — UvT°-Ul, the second term which involves the term (V^JK) gives the 
contribution I2 



h = - 



27r2 

asNc 
'(27r)2 



{x[ - X2y 



(4 - a;i)2 



{x[ — z)2(z — X2)^ {X2 — z)^{z — X\)^ 



{Tr{U{x^)U\x\)U{x2)U\x'2)\)^ 



d^z^B{xi,x[,X2,x'2;z) \ — (Tr [U{x,)U^x[)] Tr [C/(x2)C/^(4)] )^ 



1 



7V2 

^ \ ab 



{Tr[Uix^)UHx\)Uix2)UHx',)])^ 



(A2) 



^ \ ab / ^ ^ \ ab 

V^Vz) and the fourth term — ( Vj^Ki) give the contributions I3 and I4 respectively. These 



have identical forms but different variables, their sum reads 



h+h 
as Nc 
(27r)2 
a,N, 



d^z^IC2{xi;x[,x'2;z)^{TY[U{xi)U\z)]Tr[U{z)U^{x[)U{x2)UHx'2)])^ 



(27r)2 
asNc 

asNc 
'(27r)2 

asNc 



d^z^lC2{x[;xi,X2;z)^{TT[U{z)U^{x[)]TY[U{xi)U^{z)U{x2)UHx'2)])^ 



d^z^lC2{x2;x[,x'2;z)^{Tr[U{x2)U\z)]TY[U{xi)U^{x[)U{z)UHx'2)])^ 



d^z^lC2{x'2;xi,X2;z)^{TT[U{z)uHx'2)]TY[U{xi)U^{x[)U{x2)U^{z)])^ 
d'z^^{Tr[U{x,)uHx[)U{x2)UHx',)])y 



X [IC2{xi;x[,x2; z) + /C2(a;'i; xi, 2:2; z) + JC2{x2; x'^, x'2] z) + /C2(a;2; xi, a;2; z)] . 



(A3) 



Noticing that Cp — ^ ^'^'^ using Eq. it is straightforward to check that the coefhcient of the non- leading 
terms ^ (Tr [U{xi)W {x\)U{x2)W [x'^)] )^ vanishes. In the end, we find the sum of all these four contributions leads 
to the right hand side of Eq. ^ without finite Nc corrections. 
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